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INTRODUCTION 


Compared  to  the  conventional  airplane,  the  present  day  heli- 
copter is  a modern  invention.  Although  leonardo  da  Vinci  sketched  a 
man-carrying  vehicle  with  rotating  wings  around  the  year  1500,  only 
after  World  War  II  has  this  versatile  aircraft  come  into  general  use 
throughout  the  world.  During  the  last  decade  the  helicopter  has  be- 
come one  of  man's  most  useful  tools  in  the  air,  demonstrating  its  ver- 
satility in  both  civil  and  military  roles.  The  reasons  for  the  long 
delay  in  realizing  leonardo 's  dream  have  not  been  for  want  of  trying. 

If  we  consider  the  first  flight  of  de  la  Cierva's  autogyros 
in  1925-26  as  the  birth  of  rotary-wing  aviation,  then  the  "childhood 
diseases"  of  the  helicopter  were  mainly  solved  in  the  first  decade 
(1926-1936).  These  Initial  problems  were  aerodynamic  in  nature  and 
consisted  of  balancing,  controllability,  and  stability.  More  serious 
deficiencies  became  apparent,  however,  when  the  first  series-produced 
machines  appeared  and  were  placed  in  service.  These  problems  were 
dynamic  in  nature  such  as  fatigue  due  to  insufficient  dynamic  strength 
of  certain  structural  members  and  degradation  of  crew  and  passenger 
efficiency  due  to  unacceptable  vibration  levels  (Reference  1).  UTiile 
sufficient  control  of  these  deficiencies  has  been  achieved  to  allox^ 
the  helicopter  to  progress  to  its  present  state,  they  remain  as  a 
major  obstacle  in  complete  maturity  phase  development  of  the  heli- 
copter. 
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The  loads  on  the  helicopter  are  quite  different  from  a conven- 
tional propeller  driven  aircraft.  The  propeller  operates  practically 
in  an  axial  flow  and,  like  an  engine,  sets  up  no  noticeable  variable 
stresses  in  the  structural  members.  Therefore,  it  is  only  during 
take-off,  landing,  and  maneuvers  that  appreciable  dynamic  loads  are 
Induced  on  the  aircraft  structure.  Since  these  loads  usually  repre- 
sent relatively  few  load  cycles  during  the  lifetime  of  the  aircraft, 
one  can  generally  speak  of  repeated  static  loads.  The  helicopter's 
main  structural  members  are  loaded  dynamically,  the  number  of  load- 
ings often  exceeding  tens  of  millions  of  cycles  during  its  lifetime. 
The  cause  of  this  dynamic  loading  is  due  primarily  to  the  asymmetric 
airflow  past  the  main  rotor  as  it  rotates  and  simultaneously  advances 
(Reference  1) . 

In  forward  flight,  when  the  main  rotor  provides  both  lift  and 
propulsive  force  for  the  helicopter,  the  flow  over  the  blades  is  asym- 
metric due  to  the  velocity  differential  over  the  advancing  and  re- 
treating blades  (Figure  1) . Rotor  control  is  obtained  by  "cyclic 
pitch"  change,  which  is  the  name  given  to  the  first  harmonic  variation 
applied  to  the  blade  pitch  angle  as  it  rotates.  Since  the  relative 
air  velocity  over  the  blade  also  has  a first  harmonic  variation  and 
since  aerodynamic  forces  are  proportional  to  the  square  of  the  rela- 
tive velocity,  we  may  expect  to  find  at  least  three  harmonics  in  the 
force  fluctuations  acting  on  the  blades.  This  would  be  true  if  the 
airflow  through  the  rotor  were  uniform,  however,  due  to  the  proximity 
of  the  rotor  to  its  own  vortex  wake,  which  is  swept  backwards  under 
the  rotor  disc,  the  flow  is  far  from  uniform,  and  velocity  fluctua- 
tions are  Induced  which  give  rise  to  very  many  harmonics  of  blade 
loading.  The  aerodynamic  characteristics  of  a rotor  in  forward  flight 
give  rise  to  shear  forces  and  moments  at  the  blade  root  which  are  then 
transmitted  to  the  rotor  hub  where  they  are  combined  and  sent  through 
the  rotor  shaft  into  the  airframe.  The  dynamics  complexity  of  the 
helicopter  can  be  ascertained  from  the  cartoon  in  Figure  2. 

During  the  design  of  a helicopter,  accurate  prediction  of 
main  rotor  blade  and  hub  oscillatory  loading  is  important  for  fatigue, 
vibration,  and  forward  flight  performance.  Critical  dynamic  compo- 
nents are  designed  for  fatigue  based  on  these  predicted  loads.  Vi- 
bration characteristics  of  the  airframe  and  the  need  for  vibration 
reduction  devices  are  determined  from  these  loads.  Since  vibratory 
response  and  rotor  loads  usually  determine  limiting  speed  and  load 
factor  operational  envelopes,  the  predicted  loads  greatly  influence 
forward  flight  performance.  Recent  Army  helicopter  development  pro- 
grams, Utility  Tactical  Transport  Aircraft  System  (UTTAS)  and  Ad- 
vanced Attack  Helicopter  (AAH) , have  revealed  that  these  loads  can- 
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drawn  from  a comparison  study  (Reference  3)  where  the  helicopter  In- 
dustry's major  state-of-the-art  rotor  loads  analyses  were  independent- 
ly exercised  on  an  identical  hypothetical  helicopter  problem.  The 
comparison  of  results  Illustrated  significant  differences,  particular- 
ly in  structural  dynamic  modeling,  between  the  various  analyses.  A 
strong  recommendation  as  a result  of  this  study  was  to  conduct  compu- 
ter experiments  to  study  specific  Isolated  aspects  of  the  solution 
methods  and  structural  dynamics  (Reference  3) . It  was  the  results  of 
this  study  and  the  UTTAS/AAH  Programs  that  provided  the  impetus  for 
this  research. 


NIOS  lUK  AZIMUTHS  AND  VUOCnifS 


DYNAMICS  COMPLEXITY 


FIGURE  1 FIGURE  2 

THE  MATHEMATICA L MCPEL 

Description 

The  model  is  Illustrated  in  Figure  3.  It  consists  of  a 
slender  rigid  blade,  hinged  at  the  center  of  rotation,  with  spring 
restraint  at  the  hinge.  The  orientation  of  the  flap  (out-of-plane) 
and  lag  (in  plane)  restraint  springs  which  are  parallel  and  perpen- 
dicular to  the  blade  chord  line,  respectively,  simulate  the  elastic 
coupling  characteristics  of  the  actual  elastic  blade.  The  spring 
stiffnesses  are  chosen  so  that  the  uncoupled  rotating  flap  and  lag 
natural  frequencies  coincide  with  the  corresponding  first  mode  rotat- 
ing natural  frequencies  of  the  elastic  blade.  This  allows  the  model 
to  represent  a hingeless  rotor  treated  by  virtual  hinges  or  a fully 
hinged  rotor  (Reference  4) . 
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The  model  chosen  has  been  used  to  evaluate  rotor  aeroelastic 
stability  and  can  include  the  incorporation  of  flap-lag  elastic  cou- 
pling, pitch-flap  and  pitch-lag  coupling,  and  a blade  lag  damper. 
These  were  some  of  the  fundamental  dynamic  mechanisms  found  important 
in  the  stability  studies  of  References  4 and  5.  Their  Importance 
with  regards  to  forced  response  is  a major  objective  sought  in  this 
research. 

CENTRALLY  HINGED,  SPRING  RESTRAINED, 

RIGID  BLADE  REPRESENTATION 


To  Include  flap-lag  structural  coupling  and  be  able  to  vary 
it  to  determine  its  effect  on  blade  vibrations,  the  parameter  R is 
defined.  Geometrically,  it  is  defined  as  the  ratio  of  the  blade 
pitch  angle  Inboard  and  outboard  of  the  pitch  bearing.  Physically, 
it  represents  how  the  flap-lag  structural  coupling  is  dependent  on 
the  relative  stiffness  of  the  blade  segments  Inboard  and  outboard  of 
the  pitch  bearing.  This  is  because  the  principal  elastic  axes  of  the 
outboard  segment  rotate  through  an  angle  as  the  blade  pitch  varies 
while  the  inboard  segment  principal  axes  do  not.  The  method  of  im- 
plementation is  accomplished  by  replacing  6 by  R9  in  the  structural 
terms  in  the  flap-lag  equations  while  the  mass  and  inertial  terms 
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are  unchanged.  Thus,  when  R=l,  the  original  equations  are  retained, 
but  as  R is  reduced  to  zero,  the  flap-lag  structural  coupling  terms 
diminish  and  eventually  vanish  (Reference  4). 

Torsion  is  included  in  a quasi-steady  manner  as  pitch-flap 
and  pitch-lag  coupling.  This  is  accomplished  by  defining  the  pitch 
angle,  6,  as: 


0 “ Qq  ®S  + 0^6c 

Pitch-flap  coupling,  9g,  is  intentionally  built  into  some  rotors  to 
reduce  the  flapping  motion.  If  negative  6g  is  used,  then  as  the 
blade  flaps  upwardly  the  pitch  is  decreased,  which  reduces  the  aero- 
dynamic forces  on  the  blade  and  the  flapping  motion  is  somewhat  sup- 
pressed. Pitch- lag  coupling,  9^,  is  sometimes  used  in  a similar  way 
to  reduce  the  lagging  motion. 

As  stated  previously,  the  model  has  been  used  to  evaluate 
aeroelastic  stability.  These  studies,  (References  3 and  5),  have 
shown  the  model  to  give  quite  reasonable  accuracy  at  low  advance 
ratios  when  compared  with  more  complicated  models.  Since  all  pure 
helicopters  currently  developed  operate  at  ay  < 0.5  and  the  aero- 
elastic  stability  of  a rotor  is  considered  more  sensitive  to  system 
parameters  than  forced  response  the  model  is  believed  to  be  an  ade- 
quate representation  for  parametric  studies. 

Equations  of  Motion 

The  general  procedure  followed  was  to  derive  the  nonlinear 
flap- lag  equations  of  a helicopter  rotor  in  forward  flight.  The  left 
hand  side  of  the  equations  was  made  linear  by  considering  small  per- 
turbation motions  about  a periodic  equilibrium  position  (trimmed  con- 
dition) of  the  nonlinear  system.  The  nonlinear  terms  were  consoli- 
dated on  the  right  hand  side  to  form  the  following  flap-lag  perturba- 
tion equations  of  motion: 


There  are  two  forms  of  periodic  coefficients  in  the  equa- 
tions. The  explicit  periodic  coefficients  are  of  the  form  y sin  ij/ 
or  y cos  [p.  These  types  are  only  associated  with  aerodynamic  terms. 
The  second  typ^  are  Ijnplicit  periodic  coefficients  that  result  from 
the  fact  that  0 and  6 may  have  cyclic  components.  These  types  are 
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found  in  both  structural  and  aerodynamic  terms.  A damping  term  is 
included  in  the  C(ip)  matrix  to  account  for  the  lead-lag  damper  in- 
cluded on  many  rotors,  particularly  hinged  rotors.  These  equations 
are  put  into  state  variable  format  to  reduce  them  to  first  order  so 
that  an  efficient  computer  library  differential  equation  solver  can 
be  utilized. 


METHOD  OF  SOLUTION 


Eigenvalue  and  Modal  Decoupling  Analysis 

The  method  used  to  solve  the  linearized  flap-lag  equations 
with  periodic  coefficients  is  called  an  eigenvalue  and  modal  de- 
coupling analysis.  It  takes  advantage  of  Floquet  Theory  with  its  ap- 
plication to  linear  systems  with  time  varying  coefficients.  It  is 
unique  in  that  it  extends  the  Floquet  Transition  Matrix  (FTM)  Method 
(Reference  6)  to  Include  the  calculation  of  forced  response.  The 
eigenvalue  and  modal  decoupling  method  is  summarized  in  Figure  4. 

Once  the  flap- lag  equations  have  been  put  into  state  variable  format, 
the  state  transition  matrix  and  FTM  are  obtained  by  step-wise  inte- 
gration over  one  period.  This  period  consists  of  one  rotor  revolu- 
tion from  time  zero  to  2tt.  Once  the  FTM,  ♦(T),  is  obtained  it  can  be 
shovm  by  application  of  the  Floquet- llapunov  Theorem  (Reference  6) 
that  the  problem  of  stability  reduces  to  solution  of  an  eigenvalue 
problem.  Stability  is  determined  from  the  real  part  of  the  complex 
eigenvalue,  A = n + im.  Therefore,  when  n is  negative  the  system  is 
stable,  with  n = 0 representing  the  stability  limit  or  boundary. 

The  extension  of  the  FTM  Method  to  determine  forced  response 
requires  considerable  manipulation  although  the  mathematics  are  rela- 
tively simple.  The  eigenvectors  of  the  FTM  and  the  state  transition 
matrix  at  each  Integration  step  increment  must  be  obtained  and  saved. 
This  is  necessary  in  order  that  the  characteristic  functions,  A(t), 
can  be  calculated  and  used  as  a variable  substitution  to  decouple  the 
system  of  equations.  This  procedure  is  completely  analogous  to  the 
modal  analysis  method  (Reference  7)  used  to  uncouple  a system  of 
equations  with  constant  coefficients.  Hence,  the  name  eigenvalue 
and  modal  decoupling  analysis  is  given  to  the  method  of  solution. 

Hub  Shears  and  Moments 


The  information  obtained  from  the  analysis  are  the  flap  and 
lag  displacements  and  velocities.  These  quantities  are  combined  in 
the  time  dcmain  to  obtain  flap  and  lag  shears  and  moments  in  both  the 
rotating  and  non-rotating  systems.  Once  calculated  the  shears  and 
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moments  are  harmonically  analyzed  through  Fourier  Series  expansion  to 
determine  the  relative  strengths  of  the  rotor  harmonics.  In  the  ro- 
tating system  the  primary  concern  is  fatigue  of  root  end  and  hub  com- 
ponents so  that  the  first  harmonic  (1/REV)  is  the  major  oscillatory 
source.  In  the  fixed  system,  the  rotor  acts  as  a filter  and  only 
allows  rotor  harmonics  that  are  Integral  multiples  of  the  number  of 
blades  to  be  transmitted.  Therefore,  for  a two  bladed  helicopter 
only  the  even  harmonics,  2/REV,  4/REV,  etc.,  are  of  major  concern 
in  the  fixed  system.  The  root  shears  and  moments  for  each  blade  in 
the  rotating  coordinate  system  are  expressed  as  the  Integrals  of  the 
blade  aerodynamic  and  Inertial  loading.  The  single  bladed  results  in 
the  fixed  system  are  valid  for  rotors  with  any  number  of  blades  so 
long  as  the  appropriate  solution  harmonics  are  set  to  zero.  There- 
fore, if  a three  bladed  rotor  is  being  analyzed  only  the  3/REV  loads 
are  of  Interest  in  the  fixed  system. 

EIGENVAnJE  AND  MODAL  DECOUPLING  ANALYSIS 

Step  Description 

1 Put  Flap-lag  Equations  into  State  Variable  Format; 
i(t)  + D(t)X(t)  = f(t) 

2 Determine  the  State  Transition  Matrix,  <fr(t),  and 
Floquet  Transition  Matrix,  Q = i{>(T)f  By  Numerical  In- 
tegration Over  One  Period.  The  Solution  Can  Be 
Written  as:  X(t)  = (JiCt)  X(0) 

3 Find  Eigenvalues  and  Eigenvectors  of  FTM,  Q.  Floquet' s 

Theorem  States  that  a System  of  Equations  Having 
Periodic  Coefficients  Has  Transient  Solutions  of  the 
Form:  X(t)  = A(t)  At  t-0:  {a}=  A(07^X(0) 

4 Using  2 and  3 an  Eigenvalue  Problem  is  Developed: 

A(0T^Q  A(0)  = t'e^^  ].  The  Characteristic  Functions, 
A(t),  Can  Be  Obtained  From:  A(t)  = <t)(t)A(0) 

5 Through  Variable  Substitution,  X(t)  = A(t)Y(t),  the 
Equations  in  1 Can  Be  Decoupled  to  Obtain: 

^(t)  - [X]  Y(t)  = g(t) 

6 The  Uncoupled  Response,  Y(t),  Can  Be  Obtained  From 
Complex  Fourier  Series  Expansion.  The  Coupled  Re- 
sponse X(t),  Can  Be  Obtained  From  Matrix  Multiplica- 
tion. 
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7 Convert  X(t)  Information  in  Time  Domain  to  Vibratory 

Shears  and  Moments.  Fourier  Analyze  Shears 
Moments  to  Obtain  Various  Harmonic  Components. 

FIGURE  4 

RESULTS  AND  DISCUSSION 


As  with  ixi.  lew  analytical  method  correlation  with  other 
proven  analyses  or  measured  data  must  be  obtained  for  validation. 

The  hypothetical  rotor  of  Reference  3 was  used  for  this  purpose.  The 
input  parameters  for  this  articulated  (hinged)  rotor  and  for  a hinge- 
less rotor  are  presented  in  Figure  5. 


ROTOR  BUDE  INPUT  PARAMETERS 


Lag 

Type 

P 

D 

3 

^T  ^do 

R 

b c 

Ppc 

®6 

Damper 

Hinged 

1.031 

0.25  7.5  .07 

.0063  .01 

25' 

3 1.83 

0 

0 

0 

.1536 

Hinge- 

less 

1.136 

0.75  9.3  .1 

.0062  .01 

24.5’ 

4 1.92 

.06 

.15 

0 

0 

FIGURE  5 


The  first  rotor  is  typical  of  hinged  rotors  on  current  opera- 
tional helicopters.  It  includes  e lag  damper,  as  do  most  hinged 
rotors,  to  insure  freedom  from  mechanical  instability,  usually  called 
ground  resonance.  Since  the  rotor  has  three  blades  the  third  harmonic 
(3/REV)  is  the  largest  source  of  vibration  in  the  fixed  system.  How- 
ever, this  3/REV  vibration  is  caused  by  2/REV,  3/REV,  and  4/REV  shears 
and  moments  in  the  rotating  system.  The  1/REV  forces  in  the  rotating 
system  are  most  important  when  considering  blade  and  hub  component 
fatigue  lives.  The  second  rotor  is  called  hingeless  because  it  does 
not  have  a mechanical  hinge  to  allow  flapping  and  lagging  motion  as 
does  the  hinged  rotor.  The  term  "soft  inplane"  hingeless  rotor  is 
often  applied  because  the  inplane  (lag)  frequency  is  less  than  the 
rotor  speed  (1/REV).  This  rotor  includes  3pc  and  0g  in  its  baseline 
configuration.  Since  the  rotor  has  four  blades  the  first  five  rotor 
harmonics  can  be  considered  critical  in  the  rotating  system.  This 
rotor  is  very  similar  to  one  of  the  UTTAS  candidate  helicopters. 

Hinged  Rotor  Results 

Calculated  hub  vertical  and  horizontal  shears  are  compared  in 
Figures  6 and  7 with  the  results  of  Reference  3.  I'lhile  significant 
differences  exist,  the  eigenvalue  and  modal  decoupling  analysis 
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results,  labeled  FL  for  flap-lag,  compares  favorably  with  the  much 
more  sophisticated  methods. 


HUB  HOBIZONTil  SHEARS 


FIGURE  6 FIGURE  7 

Blade  lag  dampers  are  considered  essential  for  hinged  rotors 
to  Insure  freedom  from  mechanical  instability.  Coalescense  between 
the  blade  lag  motion  and  the  helicopter  rolling  or  pitching  motion  is 
the  source  of  this  instability.  Since  the  lag  motion  is  lightly 
damped  a mechanical  damper  is  added.  Figures  8 and  9 Illustrate  the 
effect  a mechanical  damper  setting  has  on  lag  damping  and  lag  moments. 
It  is  not  unusual  for  an  operating  hinged  helicopter  to  have  a 25% 
critical  damping  setting.  As  can  be  seen  from  Figure  9 this  damper 
setting  can  greatly  Increase  the  lag  moments. 

LAG  MOMENTS  » DAMPER  SETTING 
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LAG  DAMPING  vs  DAMPER 
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On  hinged  rotors  pitch-lag  and  pitch-flap  coupling  are  due  to 
mechanical  linkages  which  will  cause  the  blade  to  pitch  whenever  there 
is  a change  in  flap  or  lag  angle.  The  variation  of  flap  and  lag  damp- 
ing with  pitch-flap  coupling  is  presented  in  Figure  10.  The  effect  of 
pitch-flap  coupling  on  lag  shears  is  Illustrated  in  Figure  11.  It 
appears  that  negative  pitch-flap  coupling  can  be  used  effectively  to 
reduce  blade  loading  and  airframe  vibration  (Figure  11)  without  jeop- 
ardizing stability  (Figure  10).  The  contrary  is  true  for  positive 
pitch-flap  coupling  which  results  in  a rapid  decrease  in  flap  damping 
above  0.4. 


DAMPING  M PITCH-FUP  COUPlING  UG  SHIARS  n PnCH-FlAP  COUPlING 
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shears.  lag  damping  varies  greatly  with  pitch-lag  coupling.  The 
flattening  out  of  the  damping  curve  in  - gure  14  for  positive  pitch- 
lag  coupling  is  due  to  the  15%  critical  damping  setting  for  the  base- 
line configuration  (Figure  5). 

Hingeless  Rotor  Results 

Hingeless  rotors  are  attractive  because  they  increase  the 
helicopter's  control  power,  simplify  the  rotor  mechanically,  and  re- 
duce somewhat  the  hubweight  and  aerodynamic  drag.  With  their  hub 
cantilever  boundary  conditions  they  do  complicate  structural  dynamic 
coupling  and  transmit  large  vibratory  rotor  loads  into  the  airframe. 
Flap-lag  elastic  coupling  is  an  important  consideration  on  hingeless 
rotors  and  an  area  that  is  not  well  understood.  Figures  14  and  15 
illustrate  the  effect  of  flap-lag  elastic  coupling,  R,  on  lag  damping 
and  lag  shears,  respectively.  It  can  be  seen  in  Figure  14  that  small 
values  of  R greatly  reduce  lag  damping  but  that  the  trend  is  reversed 
for  larger  values  of  R.  With  respect  to  lag  shears;  the  1/REV  shears 
are  greatly  Increased  for  small  values  of  R and  are  relatively  con- 
stant for  larger  values. 

LAG  DAMPING  n R UG  SHEARS  vs  R 


FIGURE  14  FIGURE  15 

The  variations  with  respect  to  pitch-flap  and  pitch-lag 
coupling  are  similar  to  the  trends  found  for  the  hinged  rotor.  Figure 
16  Illustrates  the  reduction  in  flap  damping  with  positive  pitch-flap 
coupling.  For  a positive  value  of  1 the  flap  damping  is  reduced  to 
almost  zero,  the  stability  boundary.  For  positive  pitch-lag  coupling. 
Figure  17,  the  lag  damping  becomes  positive,  indicating  an  instabil- 
ity. If  no  lag  damper  had  been  installed  on  the  hinged  rotor  a 
similar  trend  would  have  been  expected. 
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FIGURE  16 
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FIGURE  17 


CONCLUSICNS  AND  RECCMMENDATimS 


An  eigenvalue  and  modal  decoupling  method  to  predict  heli- 
copter rotor  stability  and  forced  response  has  been  successfully  de- 
veloped. The  advantages  of  this  method  are: 


1.  Stability  and  forced  response  are  obtained  from  the  same  method 
thus  allowing  direct  comparison  of  parametric  effect^ 

Only  one  rotor  revolution  of  numerical  integration  for  an  initial 
condition  of  unity  imposed  on  each  degree  of  freedom  is  necessary  to 
define  the  Floquet  Transition  Matrix^  llierefore,  the  ambiguous  in- 
terpretation of  time  history  response  data  over  many  rotor  revolu- 
tions and  the  separation  of  transient  and  forced  response  is  no 
longer  necessary, 

3.  Additional  degrees  of  freedom  such  as  hub  and  airframe  motions, 
Inflow  feedback,  etc.  in  both  the  fixed  and  rotating  systems  can  be 
included  since  Floquet  theory  is  applicable  to  equations  with  peri- 
odic coefficients. 

r 

The  method  has  been  applied  to  a hinged  and  hingeless  rotor 
to  compare  the  effect  of  structural  coupling  parameters  on  stability 
and  forced  response  in  forward  flight.  Significant  conclusions  dra;,m 
from  this  application  are: 


1.  Pitch-flap  coupling  has  a powerful  influence  on  stability  and 
forced  response.  Negative  coupling  reduces  lag  shears  considerably 
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without  jeopardizing  stability,  while  positive  coupling  produces  an 
opposite  effect.  This  trend  was  similar  for  both  rotors. 

2.  Pitch-lag  coupling  has  little  influence  on  forced  response  but 
significantly  alters  lag  damping.  The  use  of  negative  coupling  re- 
sulted in  increased  lag  damping  for  both  rotors.  Positive  coupling 
reduced  lag  damping  for  the  hinged  rotor  with  lag  damper  and  produced 
an  instability  on  the  hingeless  rotor. 

3.  Variation  in  lag  damper  setting  on  the  hinged  rotor  resulted  in 
significant  increases  in  lag  moments  and,  naturally,  increased  lag 
damping. 

4.  The  elastic  coupling  parameter,  R,  altered  greatly  the  lag  damping 
on  the  hingeless  rotor.  For  low  values  of  R the  damping  decreased 
drastically  but  increased  again  for  larger  R values.  Low  values  of  R 
also  significantly  Increased  the  first  harmonic  lag  moments. 

The  current  research  has  been  directed  toward  developing  an 
alternative  solution  method  to  the  blade  stability  and  response  prob- 
lem. For  this  reason  a relatively  simple  mathematical  model  was  used 
to  study  the  effects  of  structural  coupling  parameters.  Further  re- 
search should  include; 

1.  More  case  studies  with  the  current  model  to  gain  a better  under- 
standing of  the  structural  coupling  parameters  in  forward  flight. 

Other  rotor  configurations,  such  as  "stiff  Inplane"  hingeless,  should 
be  included. 

2.  The  method  should  be  applied  to  an  elastic  blade  model  to  gain  a 
better  understanding  of  higher  harmonic  interaction. 

3.  Correlation  with  wind  tunnel  and  flight  measurements  should  be 
made  to  further  validate  the  model  and  method. 
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NOTATIONS 


a 

= 

slope  of  lift  curve 

B 

b 

= 

number  of  blades 

B 

c 

= 

blade  chord,  ft. 

Cdo 

= 

blade  profile  drag 

C(i|>) 

= 

damping  matrix 

a 

Ct 

ss 

thrust  coefficient 

6B 

KM 

B 

stiffness  matrix 

P 

= 

nondim.  rot.  flap- 

4 

ping  frequency  at 

0=0,  p=(1+u^)^ 

? 

r 

= 

blade  radial 

coordinate,  ft. 

n 

r 

nondim.  coordinate 

r/R 

0 

R 

blade  radius,  ft.. 

or  elastic  coupling 

? 

parameter 

t 

time,  sec 

T 

rotor  thrust 

®B 

Vi 

= 

uniform  induced 

velocity  in  negative 

X 

Z direction,  fps 

A 

Vx*  Vg 

components  of  heli- 

copter speed  in 

y 

negative  X and  Z 

p 

directions,  respec- 

0 

tively,  fps 

T 

X,  Y,  Z, 

aircraft  coordinates 

X,  y,  z 

and  rotating  blade 

coordinates  respec- 

tively 

X,  Y 

Forward  and  Lateral 

Shears  respectively 

L(i^) 

forcing  coefficient 

U) 

matrix 

MW 

= 

forcing  coefficient 

(03 

matrix 

KW 

SS 

forcing  coefficient 

matrix 

Q 

St(),  CiJ; 
(•) 


flap  angle,  pos  up,  rad 
equilibrium  flap  angle, 
+ BgSip  + rad 

precone  angle,  rad 
Locke  number,  pacR^/I 
perturbation  flap  and 
lag  angles,  rad 
lead-lag  angle,  pos. 
fwd,  rad 

equilibrium  lead-lag 
angle,  rad 
neg  real  portion  of 
lead-lag  eigenvalue 
pitch  angle 

e + e36B  + 036c 

equilibrium  pitch  angle 
6q  + 9gSii)  + 0 CiJ;  + 
0e(6-6pc).  rad 
pitch-flap  and  pitch- 
lag  coupling  ratios 
Inflow  ratio 
complex  eigenvalue, 

A = n+io) 

advance  ratio,  Vj^/fiR 
air  density  slugs/ft^ 
rotor  solidity,  bc/irR 
period  of  one  revolu- 
tion, T = 2iT/n  sec 
^nflow  parameter, 
t = 4/3X 

rotor  azimuth  angle, 
ip  = 0aft,  tl;  = nt,  rad, 
dimensionless  time 
imaginary  portion  of 
lead-lag  eigenvalue 
dimensionless  non- 
rotating flap  and  lag 
frequencies  at  0 = 0 
rotor  angular  velocity 
sine  (ill) , cosine  (ij<) 
_d_,  V _J . 

di>  ^ ^ * fl  dt  ^ ^ 
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